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Using firmly established experimental inputs sucli as ex, ^Md, AMs, Br(_B — > tu), 7, Vcb 
along with corresponding lattice matrix elements which have been well studied and are in 
full QCD such as Bk, SU3 breaking ratio ^, Bb,, and in particular without using Vub or the 
pseudoscalar decay constants /s^ or fn^ from the lattice, we show that the CKM-paradigm 
now appears to be in serious conflict with the data. Specifically the SM predicted value of 
sin 2/3 seems too high compared to direct experimental measured value by over 3(7. Further- 
more, our study shows that new physics predominantly effects B-mixings and Bd — >■ tpKs, and 
not primarily in kaon-mixing or in B — >■ tu. Model independent operator analysis suggests 
the scale of underlying new physics, accompanied by a BSM CP-odd phase, responsible for 
breaking of the SM is less than a few TeV, possibly as low as a few hundred GeV. Two possible 
BSM scenarios, namely warped extra-dimensions and SM with a 4th generation, are briefly 
discussed. Generic predictions of warped flavor models are briefly discussed. While SM with 
4th generation (SM4) is a very simple way to account for the observed anomalies, SM4 is 
also well motivated due to its potential role in dynamical electroweak symmetry breaking via 
condensation of heavy quarks and in baryogenesis. 

1 Introduction 

The next big step in our understanding of particle physics will be the uncovering of the electro- 
weak symmetry breaking (EWSB) mechanism. The present and upcoming collider experiments 
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(Fermilab and LHC) will be able to test the Standard Model (SM) Higgs mechanism. New 
physics is widely expected at around the TeV scale if the Higgs mass is not to receive large 
radiative corrections and require severe fine-tuning. A stringent constraint on the SM mechanism 
of EWSB is the tight structure of flavor changing (FC) interactions: tree-level FC neutral 
currents are forbidden and charged currents are controlled by the Cabibbo-Kobayashi-Maskawa 
(CKM) [T] mixing matrix 



Within the SM, the CKM matrix is the only source of FC interactions and of CP violation. 
There is no reason, in general, to expect that new physics (needed to stabilize the Higgs mass) 
at the TeV scale will be in the basis wherein the quark mass matrix is diagonal. This reasoning 
gives rise to another fundamental problem in particle physics, namely the flavor puzzle i.e. unless 
the scale of new physics is larger than 10^ TeV it causes large FCNC especially for the K — K 
system. Thus flavor physics provides constraints on models of new physics up to scales that are 
much much larger than what is accessible to direct searches at colliders such as the Tevatron or 
the LHC. Flavor physics is therefore expected to continue to provide crucial information for the 
interpretation of any physics that LHC may find. 

In the past decade significant progress was made in our understanding of flavor physics, 
thanks in large part to the spectacular performance of the two asymmetric B-factories. For 
the first time it was experimentally established that the CKM-paradigm [l| of the Standard 
Model (SM) provides a quantitative description of the observed CP violation, simultaneously in 
the B-system as well as in the K-system with a single CP-odd phase, to an accuracy of about 
20% ^|. While this success of the CKM picture is very impressive, the flip side is that an 
accuracy of 0(20%) leaves open the possibility of quite sizable new physics contributions. In 
this context it is important to recall that the indirect CP violation parameter, e/^: ~ 2 x 10^^ [s] 
is an asymmetry of O(10~^) and an important reminder that if searches had been abandoned 
even at 0(1%) the history of Particle Physics would have been completely different. Indeed, 
in the past few years as better data and better theoretical calculations became available some 
rather serious tensions have emerged (4||9]. 

Recently [lO] , we showed that the use of the latest experimental inputs along with a careful 
use of the latest lattice results leads to a rather strong case for a sizable contribution due 
to beyond the Standard Model sources of CP violation that in sin 2/3 could be around 15- 
25%. Clearly if this result stands further scrutiny it would have widespread and significant 
repercussions for experiments at the intensity as well as the high energy frontier. We also were 
able to isolate the presence of new physics primarily in the time dependent CP measurements 
via the "gold-plated" ipKs mode which intimately involves B-mixing amplitude and the decay 
B — )• "ipKs. Our analysis does not exclude possible sub-dominant effect in kaon- mixing and/or 
in B — )• TV. In particular, our analysis [10] indicates that the data does not seem to provide a 
consistent interpretation for the presence of large new physics contribution to the tree amplitude 
for B —7- TV. 

2 Choice of lattice inputs 

Key inputs from experiment and from the lattice needed for our analysis are shown in Table [T] 
Below we briefly remark on a few of the items here that deserve special mention: 




(1) 



• With regard to lattice inputs we want to emphasize that quantities used here have been 
carefully chosen and are extensively studied on the lattice for many years. Results that we use 
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A quarter century of lattice-QCD efforts to improve the determination of Bk 

(20)[2r 



As a representative of 
using lowest order ChPT and SU{3) flavor symmetry with no 



22||23 using large Nc are shown in gray. Several quenched lattice results are shown 



Figure 1: 

continuum methods, the results of Refs. 
estimate of errors and of Ref. 
in red. In particular, three earliest attempts on the lattice around 1984-85 (marked as Cabibbo et al., Gavela et 
al. and Bernard et al.) 24 -26 . Amongst these early attempts are also 27^ using Wilson Fermions and of [28| 
using staggered quarks. First large scale staggered result, marked as JLQCD is of 29 . Blum and Soni [30] marks 
the first simulation of Bk using quenched domain-wall quarks followed by large scale studies of that approach 
by CP-PACS 2001 



31 , RBC 2001 32 and CP-PACS 2008 33 



Review of all quenched results prior to 2003 is 
marked as Becirevic 34 . Unquenched 2 flavors calculations are shown in green including RBC 2004 f35 ; JLQCD 
2008 [17] and ETMC 2009 Unquenched 2-1-1 flavors calculations are shown in Blue: HPQCD 2006 [12] 

with staggered; RBC-UKQCD 2007 (38| with domain wall and ALV [T4| using mixed action. Also RBC-UKQCD 
2010 |19| us ing domain-wall quarks and SBW 16 using staggered. In black [LS] is the average used in the analysis 
of Ref. 10 and the recent average[LLV] obtained in Refs. |39lj40j. 



are obtained in full QCD with Np = 2 + 1 simulations, are in the continuum limit, are fairly 
robust and emerge from the works of more than one collaboration and in most cases by many 
collaborations. 

• Regarding calculations of I3k 
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on the lattice, it is useful to note that in the past 3 years a 
dramatic reduction in errors has been achieved and by now a number of independent calculations 
with errors < 5% and with consistent central values have been obtained using Nf = 2 + 1 (12 



16 



as well as 



17 dynamical simulations (see Ref. (18] for a review). Again, to be 



conservative, we only use weighted average of two results that are both in full QCD, use different 
fermion discretizations (domain-wall and staggered) and that also use completely different gauge 



configurations and completely different procedures for operator renormalization 16,19 



Given the important role lattice calculations of weak matrix element are playing in the 
evaluation of some of the important non-perturbative quantities that are critical to constraining 
fundamental parameters in flavor physics, we now take this opportunity to briefly comment 
on how the calculation of Bk evolved over the past ~25 years. This example should serve to 



illustrate developments in many such calculations; the history of is summarized in Fig. 
Recall that before the advent of the lattice approach to Bk , an interesting first attempt 
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21 , using charged kaon lifetime, lowest order chiral perturbation theory and flavor SU(3) sym- 



metry, estimated Bk ~ 0.33. But of course an estimate such as this represents an uncontrollable 
approximation, with no reliable error estimate or understanding of scale dependence. If one were 
to use such a value of Bk, in conjunction with experimentally measured value of ex, to deduce 
the Wolfenstein parameter rj, that uniquely controls CP violation in the CKM picture, we would 
get a central value about a factor of two higher than modern numbers, but even more noticeably 
the error on r/ could easily be 0(100%) rather than ~ 10% [T] that we now have. 

One of the primary purpose for the construction of the two asymmetric B-factories in the 
1990's was that they would allow us to extract directly from experiment the weak CP-odd phase 
via B —7- iJjKs- They accomplished this task beautifully, providing us with a rather precise 
value, sin 2/3 = 0.668 it 0.023, i.e. with an accuracy of about 3.4%. But for this accurate 
determination, attained at a significant expense and effort, to be useful in testing the Standard 
Model, and in particular to test that the CP-odd phase in the CKM paradigm is quantitatively 
responsible for the observed CP violation in Kl decays as well as in B-decays, a value for Bk with 
commensurate precision is essential. If the accuracy on Bk had stayed at the level of 0(100%) 
then the B-factory measurement would have failed to have an impact on the fundamental theory. 

The very first attempts [24-26 in the 80's on calculating Bk all started with Wilson fermions 
in the quenched approximation and of course had huge errors with a value of Bk consistent with 
or 1. 



Amongst the continuum methods, perhaps the most interesting result was that of 22 
Bk = 0.70 lb 0.10 fii], obtained by using the large N approximation: this result corresponds to 
lattice results obtained in the quenched approximation. Remarkably, this calculation has been 
consistent with all older lattice results obtained over the years in the quenched approximation, 
and in fact its claimed accuracy is higher, since, following [42] many, if not most, lattice calcula- 
tions done in the quenched approximation were quoting a (conservative) guess-estimate for the 
systematic error due to the quenched approximation of w 15% ^43j, though with hindsight we 
now see that the quenching error on Bk was less than 5%. 

In the quenched approximation, the result, Bk = 0.863 ± 0.058 (where the stated error 
does not include quenching error) of [29], obtained by using staggered quarks, and perturbative 
renormalization, stood as the most precise lattice result for a long time. With the advent of 



domain- wall quarks [30 44 45 and with large scale (quenched) simulations with domain wall 
quarks 31 ,32] it was found that domain wall quarks consistently tended to give about 8 to 15% 
smaller Bk (implying a larger CP-odd phase, rj) compared to the staggered result of [29] . 

With dynamical 2-1-1 simulations there is no longer any need for estimating quenching er- 
rors and 3-4 years ago RBC-UKQCD 13 obtained the first result in full QCD using DWQ, 
with an estimated total error of about 5.5%, finally managing to by-pass the stated accuracy 



of 23 . Furthermore, by 2010 quite a few lattice calculations using full QCD (and with dif- 
ferent discretizations flS^) have managed to reduce the error even less than that to about 4%. 
Furthermore, the latest RBC-UKQCD calculation made significant improvements in renor- 
malization and in chiral extrapolation to reduce the error further to 3.6%. 

As mentioned previously, in our study we are using a weighted average of this latest domain- 
wall result [T9I and that of Ref. [16] obtained by using staggered quarks. 



• Given the large disparity between the exclusive and inclusive determinations of at the 
level of 3.3cj, see Table [H it is very difficult to draw reliable conclusions by using this quantity; 
therefore, since 2008 |5| we have been advocating not using V^b for testing or constraining the 
UT. Consequently in this work also we will make very limited and peripheral use of Vub only. 
We should also stress that this is one of the key differences between other groups [61 [9] work 
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lK6lexci = (39.5 ±1.0) X 10-^ 
iV^ftli^^i = (41.68 ± 0.44 ± 0.09 ± 0.58) x lO'^ |50 
l^cfelavg = (40.9 ±1.0) X 10-3 
|Kfelexci = (31-2 ±2.6) X 10-4 
|Kfeli,,i = (43.4 ±1.611;^) X 10- 
iKftltot = (33.7 ±4.9)10-4 
AruBa = (0-507 ± 0.005) ps-^ 
AruB, = (17.77 ±0.12) ps-^ 
EK = (2.229 ±0.012)10-3 
mt,poie = (172.4 ±1.2) GeV 
rricimc) = (1.268 ± 0.009) GeV 
S^Ks = 0.668 ± 0.023 [Ssj 

/b,/^= (291±16) MeV 
BRb^^^ = (1.68 ± 0.31) X 10-4 



61-63 



r]i = 1.51 ± 0.24 49 
rj2 = 0.5765 ± 0.0065 
r/3 = 0.494 ± 0.046 
7]B = 0.551 ±0.007 
^ = 1.23 ±0.04 
A = 0.2253 ± 0.0009 
a = (89.5 ±4.3)° 
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Ke = 0.94 ± 0.02 
Bd = 1.26 ±0.11 
fB, = (208 ± 8) MeV _ 
fx = (155.8 ± 1.7) MeV 
Bk = 0.742 ± 0.023 
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7 = (78± 12)° 
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Table 1: Lattice QCD and other inputs to the unitarity triangle analysis. The determination of a is obtained 
from a combined isospin analysis of B 



(tttt, pp, pTv) branching ratios and CP asymmetries 50 . Statistical and 
systematic errors are combined in quadrature; for the error on Vub see [64]. We adopt the averages of Ref. |39||40| 
(updates at www.latticeaverages.org) for all quantities with the exception of ^, /b^bI^^, Bk and fs^ (see 
text) . 



^1/2 '^1 /2 

on UT fits and us. Another difference from those works is that we use ^ = /bs^s / fsd^d ' 
IbsBs and B^ to describe Bg mixing and B ^ tv (as opposed to fBsl fB^-, Bg/B^, Bg and 
/Bs)- Moreover, we fit /b^ in conjunction with particular hypotheses and use the determined 
value of as a diagnostic tool. Another difference between our work and Ref. j9j is that, in 
the latter, the authors include in the fit the DO dimuon asymmetry 
Bg and Bg decays; while we agree that this is a very interesting result, we believe that it needs 
confirmation and therefore are not using it in our fit. 
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47| , agj^ , in semileptonic 



• The complete set of lattice inputs that we use is presented in Table [T] All inputs, are taken 
from Refs. [39y40| (see http://www.latticeaverages.org for updates) with the exception of 
Bk (see discussion above), ^ (since the statistical errors of the HPQCD and Fermilab/MILC 
results are 100% correlated, we decided to increase the statistical error of the HPQCD result to 
bring it in line with the with the more conservative Fermilab/MILC estimate), /b^ (we update 

; " ^1/2 

the HPQCD determination of /^^ 148]) and /bs^s (we update the HPQCD determination of 



/Bs 48 and combine it with the Fermilab/MILC result; we then combine the /b^ average with 



the HPQCD determination of Bg adding linearly the uncertainties.) 



3 Some results of the fit. 

We first draw attention to the results of the fit shown in the upper panel of Fig. [2] Here we 
use as inputs from experiments, ex, ^M^, AM^, 7 and BR(i? — tu) [65j and from the lattice, 
i^K, fBsBl/2 and I3ii (but not /b^) and we extract the fitted value of sin 2/3 and of Jb^- We 
obtain: 

sm{2(3f^ = 0.867 ± 0.050 , (2) 



which is about 3.2 a away from the experimentally measured value of 0.668 ± 0.023. We believe 
this result provides a strong indication that the CKM description of the observed CP violation 
is breaking down 66 . 



For the fitted value of /b^ along with the predicted value of sin(2/3) given above, we find: 



f^l = (201.5 ± 9.4) MeV . (3) 

This "predicted" value is in very good agreement with the one obtained by direct lattice calcula- 
tion, /b^ = (208 ± 8) MeV. This is a useful consistency check signifying that the SM description 
of the inputs used, especially of i? — )• rz^, is working fairly well and that it is unlikely that 
the B ^ Tu tree amplitude is receiving large contributions from new physics; most likely the 
dominant effect of new physics is in fact in sin(2/3). Later we will reexamine this from an en- 
tirely different perspective and show in fact there is additional independent support to these 
interpretations. 

In order to further scrutinize the tentative conclusion reached above, we next present an 
alternate scenario depicted in the bottom panel of Fig. [2| Here, we make one important change 
in the inputs used. Instead of using the measured value of BR(i? — >• tv) we now use as input 
the measured value of sin(2/3) from the "gold-plated" i^Kg mode. Again, this fit yields 

two important predictions: 

BR(5 ^ Tvf^ = (0.768 ± 0.099) x 10"^ , (4) 
/i*= (185.6 ±9.1) MeV. (5) 

Eq. Q deviates by 2.7a" from the experimental measurement, as can also be gleaned from an 
inspection of the bottom panel of Fig. [2| It is particularly interesting that also the fit prediction 
for fB^ now deviates by about I.80" from the direct lattice determination given in Table [ij We 
believe this provides additional support that the measured value of sin(2/3) being used here as a 
key input is not consistent with the SM and in fact is receiving appreciable contributions from 
new physics. 

This leads us to conclude that while the presence of some sub-dominant contribution of new 
physics m. B ^ TV is possible, a large contribution of new physics in there is not able to explain, 
in a consistent fashion, the tension we are observing in the unitarity triangle fit. 

This conclusion receives corroboration by the observation that even without using B ^ tu 
at all, and using as input only ek, ^MbJ^Mb^ and \Vcb\ (see Fig. [s]), the predicted value of 
sin(2/3) deviates by 1.8cj from its measurement (in this case we find sin(2/3)'^* = 0.814 ± 0.081). 
Thus, possible new physics va. B ^ tu can alleviate but not remove completely the tension in 
the fit. 

We recall that the fit above is actually the simple fit we had reported some time ago (now 
with updated lattice inputs) with its resulting ~ 2 o" deviation [sj. This fit is somewhat special 
as primarily one is only using AF = 2 box graphs from ek and AMs^/AM^^ in conjunction 
with lattice inputs for Bk and the SU(3) breaking ratio ^. The experimental input from box 
graphs is clearly short-distance dominated and for the lattice these two inputs are particularly 
simple to calculate as the relevant 4-quark operators have no mixing with lower dimensional 
operators and also require no momentum injection. The prospects for further improvements in 
these calculations are high and the method should continue to provide an accurate and clean 
"prediction" for sin(2/3) in the SM. So even if the current tensions get resolved, this type of 
fit should remain a viable way to test the SM as lattice calculations and experimental inputs 
continue to improve. 

3.1 Roles ofVcb, £k, Vub ctnd of hadronic uncertainties. 

The fit described above does use Vcb where again the inclusive and exclusive methods differ 
mildly (about 1.7cj). Of greater concern here is that ex scales as |V^fe|^ and therefore is very 
sensitive to the error on Vcb. We address this in two ways. First in the upper panel of Fig. [3] we 
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Figure 2: Unitarity triangle fit. In each plot inputs that are grayed out are not used to obtain the black contour 
(which represents the SM allowed la region), the p-value and the fit predictions presented in the upper left 
corners. The deviations of the fit predictions for sin(2/3) and BR(i? — >■ tv) from the respective measurements 
are obtained using the actual chi-square distribution for these quantities. The p-value of the complete SM fit 
(i.e. including all the inputs) is psM = 1.7%. In the upper panel, we consider a scenario with a new phase in Bd 
mixing, thereby removing the sin(2/3) and a inputs. In the lower panel we consider a scenario with new physics 
in B — >■ TV, thereby removing the BR(i3 — )■ tv) input. 



study a fit wherein no semi-leptonic input from 6 — t- c or 6 — t- n is being used. Instead, in this 
fit BR(i? — )• Tu) and AM^^ along with ^Mb^/AMb^ and 7 are used. Interestingly this fit 
gives 

sin(2;g)^* = 0.905 ± 0.047 , (6) 
/I* = (202.9 ± 9.3) MeV (7) 

Thus, once again, sin(2/3) is off by 3.1a whereas /b^ is in very good agreement with directly 
measured value which we again take to mean that the bulk of the discrepancy is in sin(2/3) rather 
than in B ^ TV or in Vcb- 

Next we investigate the role of ck- In the bottom panel of Fig. [3] we show a fit where only 
input from B-physics, namely AMb^/AM^^, AMb^, 7, |Vc6| and BR(S — )■ ri^) are used. This 
fit yields, 

sin(2/3)*^* = 0.889 ± 0.055 , (8) 
/fit = (200.7 ± 11) MeV (9) 

Thus, sin(2/3)fi* is off by w 2.4cr and again is in good agreement with its direct determination. 
We are, therefore, led to conclude that the role of ex in the discrepancy is subdominant and that 
the bulk of the new physics contribution is likely to be in B-physics. As before, the fact that 
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Figure 3: Unitarity triangle fit without semileptonic decays (upper panel) and without use of K mixing (lower 
panel) . See the caption in Fig. |2] 

the fitted value of fs^ is in good agreement with its direct determination seems to suggest that 
the input BR(i3 — t- tu) is most hkely not in any large conflict with the SM, though, obviously 
we cannot rule out the possibility of it receiving a sub-dominant contribution from new physics. 

For completeness, we present in Fig. [4] the results we obtain when including V^b in the fit. 
Note that inclusive and exclusive determinations of \Vuh\ differ at 3.3(T (see Table[T]) and, for this 
reason, are presented separately in the plot. Before taking the average, we add a 10% model 
uncertainty to the inclusive determination. This reduces the discrepancy to 2.1a". We finally 
rescale the error on the average by the square root of the reduced chi-square (following the PDG 
recipe). In Table [l] we report the result we obtain and that we use in the fit. 

A compilation of all the eleven fits that we studied for sin 2/3 are shown in Fig. [5| Notice that 
there is only one case in here (8th from the top) where the discrepancy in sin 2/3 is only 0{la). 
We believe this is primarily a refiection of the large (~ 14.4%) uncertainty with our combined 
Vub fit originating from the large disparity between inclusive and exclusive determinations. This 
is again a reminder of the fact that till this discrepancy gets removed, we cannot use Vub to draw 
any reliable conclusion. 

4 B ^ Tu and new physics 

Now with regard to S — )• tv, Fig. [6] shows a summary of predictions versus the measured BR. 
Notice that whenever the measured value of sin(2/3) is used as an input, the predicted BR is 
~ 2.7(7 from the measured one. In the preceding discussion we have emphasized that this seems 
to us to be a consequence of new physics largely in B mixings and/or in B^ — t- ^IjKs decay. This 




conclusion receives further strong support when we try determine the B ^ tv branching ratio 
without using sin 2/3. Indeed as shown in Fig. |6] when we use eK, AM^^, Vcb and 7 only, the 
fitted value of BR(i3 — t- tu) is in very good agreement with the measured value. 

In principle, of course, the prediction for BR(i? — t- tv) only needs the values of Jb^ and 
of Vub- Fixing now fs^ = 208 ± 8 MeV as directly determined on the lattice (see Table [T]) we 
show the corresponding two predictions for the BR using separately the values of V^b determined 
in inclusive and in exclusive decays. It is clear that the inclusive determination yields results 
that are within one a of experiment (see also Fig. [2]); however with V^b from exclusive modes 
(that makes use of the semileptonic form factor as determined on the lattice), the BR deviates 
by ~ 2.8(T from experiment. This may be a hint that lattice based exclusive methods have 
some intrinsic difficulty or that the exclusive modes are sensitive to some new physics that 



the inclusive modes are insensitive to, e.g. right-handed currents 67,68 . In either case, this 



reasoning suggests that we try using the value of Vub given by inclusive methods only in our fit for 
determining sin 2/3. This line of reasoning is also supported by the analysis presented in Ref. [69| 
in which the discrepancy between the experimental determination and the SM prediction of ratio 
Rgji = BR(i? — )• niu) /BR(i? — )• tu) is considered. Note that the authors of Ref. 69 find that 
the experimental value of this ratio is about a factor of 2 smaller than the SM prediction and 
that this discrepancy is independent of whether lattice QCD or Light-Cone QCD Sum Rules are 
used to determine the i? — )• vr form factor and the B decay constant. This result can be seen as a 
solid consistency check of the lattice QCD calculation of the B ^ it form factor. Within the SM 
this ratio is independent of short distance physics (the \Vub\'^ factors cancel out) and measures 
the ratio of the -B — ?■ vr form factor to the B decay constant. New physics in right-handed 
currents affects differently the B — t- Triv and B ^ tv transitions and might be responsible for 
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Figure 5: Summary of sin(2/3) determinations. The entry marked *** (tenth from the top) is obtained with 
lattice errors increased by 50% over those given in Table [l] for each of the input quantities that we use and the 
entry marked +-1—1- (eleventh from the top) corresponds to adding an hadronic uncertainty &^S^k ~ 0.021 to 
the relation between sin(2/3) and S^k- See the text for further explanations. 
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the observed discrepancy. 



5 Summary of fits, perspective &: outlook. 

The result of our analysis strongly suggests that the SM predicted value of sin(2/3) is around 
0.85 whereas the value measured experimentally via the gold plated ^l^Ks mode is around 0.66 
constituting a deviation of about 3o" from the SM (see Fig. [s]) . To put this result in a broader 
perspective let us now recall that in fact in the SM sin(2/3) can also be measured via the 
penguin dominated modes (see Fig. [s]) 70 - 73 . Unfortunately several of these modes suffer from 



a potentially large tree pollution, though there are good reasons to believe that the r]'Ks, (jjKs 



and 3 Kg modes are rather clean 74-76 wherein the deviations from sin 2/3 are expected to be 
only 0(few %). The striking aspect of these three clean modes as well as many others penguin 
dominated modes (see Fig. [s]) is that the central values of almost all of them tend to be even 
smaller than the value (0.66), measured in ipKs, and consequently tend to exhibit even a larger 
deviation from the SM prediction of around 0.85. Thus, seen in the light of our analysis, the 
deviation in these penguin modes suggests the presence of new CP-violating physics not just in 
S-mixing but also in 6 — )• s penguin transitions. 

Moreover, the large difference (~ (14.4 it 2.9)%) js] in the direct CP asymmetry measured in 
—7- K^-K~ versus that in — )• K'^'k^ provides another hint that 6 — )• s penguin transitions 
may be receiving the contribution from a beyond the SM source of CP-violation (for alternate 



explanation see Refs. 77-79]). To briefly recapitulate, in the SM one naively expects this 
difference to be vanishingly small and careful estimates based on QCD factorization ideas suggest 
that it is very difficult to get a difference much larger than (2.2 it 2.4)% [Tj. 

Of course, if 6 — )• s penguin transitions (AFlavor = 1) are receiving contributions from new 
physics, then it is quite unnatural for Bs mixing amplitudes (AFlavor = 2) to remain unaffected. 
Therefore, this reasoning suggests that we should expect non-vanishing CP asymmetries in 
Bg — )• ip(f) as well as a non- vanishing di-lepton asymmetry in Bs — ?• Xglv. As is well known, 
at Fermilab, in the past couple of years CDF and DO experiments have been studying CP 
asymmetry in Bs — )• V"^- The latest result with about 6 fb~^ from each experiment seems 



to reveal a reduction from ~1.8cr tension to ~lcr from the SM 50,80,81 . Thus, findings in 
Bs — 7- V"^ from Fermilab and from LHCb are eagerly awaited. 

Another interesting and potentially very important development with regard to non-standard 
CP in Bs is that last year DO announced the observation of a large dimuon asymmetry in B- 
decays amounting to a deviation of (~ 3.2(t) from the minuscule asymmetry predicted in the 
SM [46|47| . They attribute this largely to originate from Bs mixing. While this is a very exciting 
development, their experimental analysis is extremely challenging and a confirmation is highly 
desirable before their findings can be safely assumed. Note, though, HFAG [50j has combined 
CDF and DO results on Bs — )• ip4> and on the dimuon asymmetry, and finds the deviation 
from the SM to be around 2.7cj. 

Be that as it may, we reiterate that our analysis suggests that the deviation from the SM in 
sin(2/3) is difficult to reconcile with errors in the inputs from the lattice that we use, and strongly 
suggests the presence of a non-standard source of CP violation largely in B/Bs mixings, thereby 
predicting that non-standard signals of CP violation in S{Bii — )• r]'Ks,(f)Ks,3Ksetc.) as well as 
in S{Bs — )• V"^)) Slid the semileptonic and di-lepton asymmetries in Bs, and possibly also in B^, 
decays will persist and survive further scrutiny in experiments at the intensity frontiers such as 
Fermilab (CDF, DO), LHCb and the Super-B factories. Lastly, the fact that our analysis rules 
out the possibility that new physics exclusively in kaon mixing is responsible for the deviations 
in sin(2/3), has the very important repercussions for the mass scale of the underlying new physics 
contributing to these deviations: model independent analysis then imply that the relevant mass 
scale of the new physics is necessarily relatively low, i.e. below 0(2 TeV) [7,82 . Thus, collider 



experiments at the high energy frontier at LHC and possibly even at Fermilab should see direct 
signals of the underlying degrees of freedom appearing in any relevant beyond the Standard 
Model scenario. 



6 Aftermath: BSM possibilities 

Let us next discuss a model independent point of view as to the possible underlying cause for 
these anomalies and then two specific models that may be relevant. 



6. 1 Brief Summary of the model independent analysis 

One of the important issue is how these B-CP anomalies will impact search for New Physics at 
the LHC wherein a knowledge of underlying scale of NP would be very useful. With this in mind 
we [T] write down dimension-6 operators under the general assumptions of NP in AFlavor = 2 
effective Hamiltonian for K, or Bg mixing or for the case of AFlavor = 1 Hamiltonian that 
may be relevant for 6—7-5 penguin transitions |7]. Our model independent analysis shows that 
the scale of CP violating NP is only a few hundred GeV if it originates from 5—7-5, AFlavor = 1 
penguin Hamiltonian. It rises to about a few TeV if it originates from B^ and/or Bg mixing. 
From the perspective of LHC the scenario that is the most pessimistic, with NP scale in the 
range of a few tens of TeVs, is when all of the NP resides only in the dimension-6 LR-operator 
relevant for the K — K mixing j83| . However, in the preceding sections we have shown that bulk 
of the deviation from the SM does not originate va. ex ox K — K mixing. This, therefore, has 
the important consequence that taken seriously these discrepancies with the CKM hint at scale 
of new physics that is quite likely to be less than a few TeV, possibly even a few hundred GeV. 
Note also that in our 08-09 work |7) we were unable to rule out the possibility that dominantly 
NP resides in /C-mixing; this only became possible in our 2010 analysis |10|. 



6.2 Warped Flavordynamics & duality 

Perhaps the most interesting and even compelling BSM scenario is that of warped extra dimen- 
sional models |84| as it offers a simultaneous resolution to EW-Planck hierarchy as well as flavor 



puzzle [85 86 . While explicit flavor models are still evolving, potentially this class of models has 
many interesting features; for example, in general one expects several new BSM CP-odd phases 
(presumably 0(1)) [8^ that can have important repercussions for flavor physics. Indeed in the 
simplest scenario it was predicted [88] that there should be smallish {i.e. 0(20%)) deviations 
from the SM in B^ decays to penguin-dominated final states such as (pKg, ij'Kg etc as well as 
the possibility of a largish CP-odd phase in Bg mixing which then of course has manifestations 



in e.g. Bg — ?- time-dependent asymmetries in Bd — )• {p,K*)j, etc ^,88 . In fact there was 



also a mild CP problem in that very simple rough estimates suggest neutron-EDM should be 



bigger than current bounds 89 by about an order of magnitude. 

Note though that in this original study, for simplicity, it was assumed that B^, mixing was 
essentially described by the SM. More recently there have been two extensive studies of the 
possibility of warped models being the origin of the several hints in B, Bg decays mentioned 
above [9^|92]. 

A common feature of these warped models is that they also imply the existence of various 
Kaluza-Klein states, excited counterparts of the gluon, weak gauge bosons and of the graviton 
with masses heavier than about 3 TeV |93|. Note also that unless the masses of these particles 



are less than about 3 TeV their direct detection at LHC will be very difficult [94 -102|. 



Since viable explicit models are still being developed, it is useful to emphasize generic pre- 
dictions of warped flavor models. Top quark should exhibit large flavor violations via e.g. 



t —7- c{u)Z 102 , largish mixing with possibly enhanced CP violations, beyond the SM CP 



asymmetries are also possible in B^, Bg physics, also polarized top quarks endowed with some 
forward-backward asymmetry |103| should be expected. Furthermore, KK particles such as the 
KK-glue, graviton, Z' should have large BR to top quark pairs and in fact the tops should be 



boosted since the decaying resonances are expected to have TeV-scale masses 94 - 97 , 99 



An extremely interesting subtlety about these 5-dimensional warped models is that they are 

This serves 



supposed to be dual to some 4-dimensional models with strong dynamics 104 -110 



as motivation to search for effective 4-dimensional models that provide a good description of the 
data. 



6.3 Extension of SM to four generations: SM4 

SM with four generations provides a rather simple and an interesting extension to address the 



B-CP anomalies [111-118 . Actually,for several reasons, SM4 is of considerable interest irrespec- 



tive of these anomalies: 

• The heavier quarks could form condensates and thereby play an important role in dynamical 
electroweak symmetry breaking [119 -122 . 



Two new CP -odd phases and the heavier quark masses also significantly alleviate one of the 



key difficulty that SMS has for baryogenesis 123 



• SM4 can open new avenues for addressing the dark matter issue 1124 125 



• Besides, since we have already seen three families, it is natural to ask why not the fourth? 

Note also that while practically all BSM scenarios have difficult time accounting for the 
(almost) absence of FCNC (in processes such as BR(i? — )• Xg^)), SM4 explains this readily. First 



of all in SM4 (as in SM) FCNC are loop suppressed via the GIM [126 mechanism. Furthermore, 



as you extend the 3x3 matrix to 4 x 4 and impose unitarity, the hierarchical structure of CKM 
matrix extended to 4 x 4 automatically allows only small residual CP-conserving effects in 
quantities such as BR(i? — )• Xg'y); on the other hand, there can be dramatic difference in CP 
violating observables where in the SM null results are predicted [127 . 



In contrast to CP-conserving FCNC, since there are now two new CP-odd phases, they can 
cause large (0(1)) CP -asymmetries in channels that the CKM phase has negligible effect in 
the SM [127j . This is expected to be the case e.g. in Bs mixing, causing S{Bs — )• ipcp) and 
the semileptonic asymmetry, a*j, in Bg — )• Xgii' to be non- vanishing. Similarly there may be 
non-vanishing CP-asymmetries in 6 — t- 57, b — t- sll, B^ — t- rj'Kg, B^ — s- (pKg wherein the SM 
one expects negligible effects. Moreover, we should also expect the effects due to an additional 
phase in B^ mixing beyond what's there in the SM. This manifests itself say as a deviation in 
e.g sin 2/? from the SM predicted value and also can cause the semi-leptonic asymmetry, a^^ to 
differ from the SM predicted value (that is negligibly small). 

It may be useful also to note that seen from the perspectives of SM4 the hierarchy puzzle 
may just be a historical accident. 

Natural mass scale for a Higgs particle in SM4, where the heavy quarks are geared towards 
EW symmetry breaking, is around 2m;,/ . Such a heavy Higgs would of course have very clean 
decays to H ^ ZZ. 

Electroweak precision tests do not rule out the existence of a 4th generation though they 
restrict the mass splitting between the 4th generation doublet {t',b') of quarks to be less than 
around 75 GeV [128 , 129 . This requires some 10% degeneracy in their masses. One cautionary 
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Figure 7: Correlation between SipKs and the semi-leptonic asymmetry, af; in SM4 is shown for mfi varying 
between 375 and 575 GeV. The experimentally allowed region as well as the SMS bounds are also shown. See 
Ref. [mj for details. 



remark is that such heavy quark masses means large Yukawa coupHngs, therefore many pertur- 
bative calculations, including those relevant for EW precision tests may receive large corrections. 

Furthermore LEP experiments require that the 4th generation neutral lepton has to be rather 
heavy > mz/2; this begs the question as to why there should be such a huge disparity with 



the three almost massless neutrinos of the conventional SMS 124 , 125 . Issues such as these are 
interesting and require further investigations. 

We will briefly now summarize the possible relevance of SM4 to alleviating the discrepancies 
in the CKM picture that has been our primary focus here. In particular, the tension in S{^l^Ks 



that we emphasized has been examined in 113 . Fig from Ref. 113 shows a study of the 



correlation between the predicted value of sin 2/3 and a^^ , which is the semi-leptonic asymmetry 
in Bel ~^ Xdlu. In this study all known experimental constraints such as ex, Br{K^ — t- vr+i/p), 
AMf^, AMs, constraints on the unitarity angle 7 etc have all been incorporated. As is evident 
from the figure, while the SM4 can accommodate the measured value of '&{ipKs) it also requires 
simultaneously that a^^ > —0.001, which is only about a factor of a few different from the SMS. 
This underscores another attractive aspect of SM4: it is rather predictive, highly constrained 
extension of the SM and can be ruled out with relative ease. In this specific illustration it 
requires improved determination of a^^ as well as S^k^ both theoretically and experimentally. 
Improved bound from the B factories, who need to update their several years old result [50] on 
a^; would be very useful. 

Another interesting example is the semi-leptonic asymmetry in Eg, a'^i and its correlation 
with S{Bs — )• V"/*)) see Fig. [sjfrom 113 . It is interesting to note here that while SM4 allows a 



much larger semi-leptonic asymmetry as shown, it is still not large enough to explain the central 
value of the asymmetry reported by the recent DO result |47) . Thus, if improved experimental 
results uphold near the central value of the DO experiment, then SM4 may also be ruled out. 

Finally, let us also briefly mention that experimental searches for quarks (t', b') of the 4th 
generation have already been underway at Fermilab leading to a lower bound of around 350 
GeV [ISO . We should expect intensified searches at LHC especially since this is something that 



can be achieved even in the early 7 GeV run 131 . It is expected that after several years of 



efforts, LHC should be able to find these quarks or put a bound close to a TeV 132 
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Figure 8: Correlation between SijjKs and the semi-leptonic asymmetry, af; in SM4 is shown for varying 




between 375 and 575 GeV. The experimentally allowed region (la) is shown. See Ref. 113 for details 
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